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HEAT  TRANSFER  STUDIES  IN  THE  REGION  OF  SHOCK  AND  TURBULENT 
BOUNDARY  LAYER  INTERACTION  INDUCED  BY  A  CYLINDRICAL  PROTUBERANCE* 


< 


Yu  Hong-ru,  Li  Zhong-fa 

|  (Institute  of  Mechanics  Academia  Sinica) 

A 

t  SUMMARY 

■j  This  article  describes  the  results  of  experiments  on  the 

rates  of  heat  exchange  in  areas  of  interference  associated  with 
concave  cylindrical  objects.  In  these  experiments,  it  has 
appeared  that  the  maximum  rates  of  heat  transfer  on  the  center 
line  of  the  upper  flow  do  not  occur  at  the  foot  of  the  cylindri¬ 
cal  object  itself  but  rather,  at  a  place  about  0.1  to  0.15  times 
the  length  of  the  diameter  of  the  cylinder  from  the  base  of  the 
cylinder.  In  the  case  of  a  vertical  cylinder,  when  M  =  9,  the 
largest  rates  of  heat  flow  are  46  times  the  corresponding  values 
for  the  same  area  when  a  flat  plate  Is  used  Instead  of  the  cylinder. 
Moreover,  when  M  *  6.6,  these  values  are  31  times  the  correspond¬ 
ing  values. 

The  cylinder  inclination  has  a  very  large  influence  on  the 
Interference  areas  produced  and  on  the  maximum  rates  of  heat  trans¬ 
fer.  This  becomes  most  obvious  when  the  Mach  numbers  Involved  are 
high.  When  one  Is  dealing  with  a  backward  sweep  of  30°,  the  maxi¬ 
mum  rate  of  heat  transfer  is  1/6.  A  forward  sweep  causes  an 
increase  in  both  the  area  of  interference  and  the  maximum  rate  of 
heat  transfer. 

1 .  FOREWORD 


The  surfaces  of  high  speed  craft  and  aircraft  cannot  avoid  the 
necessity  for  the  installation  of  various  types  of  measuring  devices, 
antennas  and  other  such  similar  protuberances  as  tail  fins.  The 


This  article  was  received  on  January  14,  1980. 


1 


blunt  forward  edges  of  these  protuberances  produce  interference 
between  the  surface  boundary  layers  of  the  aircraft  or  spacecraft 
involved  and  the  bow-shaped  detached  shock  waves.  The  occurrence 
of  flow  separation  and  reattachment  causes  a  great  increase  in 
the  heating  rate  in  the  neighboring  areas  as  compared  with  a  sit¬ 
uation  in  which  there  is  no  protuberance  and  this  is  particularly 
true  when  the  Mach  numbers  involved  are  high.  Due  to  the  com¬ 
plexity  of  this  type  of  three-dimensional  shock  wave-boundary  layer 
interference  flow  field,  it  is  first  necessary  to  make  experimental 
observations  and  measurements  in  order  to  make  clear  the  principles 
which  govern  these  interference  flow  fields  and  their  new  pheno¬ 
mena.  This  sort  of  initial  work  makes  it  easier  to  construct 
appropriate  flow  models. 

Concerning  the  typical  sorts  of  cylindrical  protuberences ,  this 
matter  has  drawn  the  attention  of  quite  a  few  researchers.  However, 
most  of  them  have  concentrated  on  the  clear  representation  of  the 
surface  flow  spectra  caused  by  turbulence  as  well  as  on  the  area 
influenced  by  this  turbulence  [1-5].  Turbulence  flow  interference 
heat  transfer  experiments  have  been  limited  to  the  supersonic  realm 
[6].  The  Mach  number  of  the  gas  flow  involved  has  a  very  important 
influence  on  the  interference  area  heat  transfer.  High  Mach  number 
three-dimensional  shock  wave  turbulence  flow  boundary  layer  inter¬ 
ference  heat  transfer  experiments  are  very  important. 

The  area  of  the  most  severe  Increase  in  heat  is  located  in  the 
base  of  the  cylinder.  The  heat  flow  distribution  at  this  point 
changes  dramatically.  If  one  does  not  carefully  position  the  point 
at  which  measurements  are  taken,  it  is  difficult  to  measure  accur¬ 
ately  the  peak  values  and  the  positions  of  these  values.  This 
article  puts  its  main  emphasis  on  the  research  into  the  heat  trans¬ 
fer  distribution  characteristics  of  this  area  of  severe  increases 
in  heat . 
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In  the  winter  of  1977,  when  we  were  running  experiments  at 
Mach  number  M  =  9  conditions,  we  made  complete  measurements  of 
the  heat  flow  rates  on  the  surfaces  around  cylindrical  protuber- 
ences.  We  found  that  the  maximum  surface  heat  flow  rates  on  the 
center  line  of  the  upper  flow  of  the  cylinders  considered  were 
located  at  a  point  approximately  1/10  the  diameter  of  the  cylinder 
from  the  surface  of  the  cylinder.  The  maximum  values  which  were 
measured  were  far  higher  than  the  numerical  values  which  were 
obtained  by  calculations  based  on  empirical  formulae  [7,8]  cur¬ 
rently  in  use  in  China  and  overseas.  Due  to  the  fact  that,  at 
that  time,  we  were  using  plug  type  sensors,  these  sensors  suffered 
from  the  limitations  of  their  dimensions  (diameter  5.5mm);  more-  ! 

over,  the  placement  of  these  sensors  was  relatively  sparse.  These 
experiments  made  use  of  changes  in  the  diameter  of  the  cylinders 
involved  and  results  were  only  obtained  after  non-dimensional  pro¬ 
cessing  with  the  diameter  of  the  cylinder  involved  as  the  basis. 

In  order  to  verify  the  accuracy  of  this  set  of  principles  for 
describing  the  heat  transfer  parameter  distributions  involved  as 
well  as  the  accuracy  of  the  data  offered,  a  new  series  of  tests 
was  run  In  1978  after  the  test  production  of  a  new  sensor  (width 
0.5mm,  sensor  Interval  2mm).  Due  to  the  relative  density  of  the 

i 

distributed  sensor  points,  the  maximum  values  which  were  measured 
for  the  surface  heat  flow  rates  were  approximately  20%  higher  than 
those  measured  during  the  previous  tests.  The  maximum  surface  heat  ; 

flow  rates,  compared  to  those  measured  when  there  was  no  concave  5 

protuburence  present  and  measured  for  the  same  positions  with  a  • 

flat  plate  Installed,  were  46  times  higher.  In  order  to  observe 
the  influence  of  the  gas  flow  Macb  number,  in  1979,  another  set  of 
supplementary  experiments  were  carried  out  with  M  «  6.6.  This 
article  gives  an  overall  description  of  the  results  of  these  exper¬ 
iments  with  special  emphasis  on  the  various  types  of  phenomena  which 
appeared  in  these  experiments. 


2.  A  DESCRIPTION  OF  THE  EXPERIMENTS 


Experimental  equipment.  These  experiments  were  carried  out 
in  the  JF8  shock  wave  wind  tunnel.  The  test  section  was  6.5m 
long,  and  the  test  length  was  13.3m.  The  internal  diameter  in 

both  cases  was  150mm.  The  semi-angle  for  the  cone-shaped  jet  tube  !  » 

1  I  * 

was  7— °.-  •  The  jet  tube  exhaust  diameter  and  the  diameter  of  the  i 

8  i 

experimental  section  were  the  same,  1.2m. 

The  experimental  gases  were  air  and  use  was  made  of  two  types 
of  operational  conditions.  The  first  of  these  was  a  reverse  jet 
model  of  operation.  The  total  pressure  in  the  storage  chamber  was  .  ;  j 

420  bars.  The  total  temperature  was  950°K,  and  the  Mach  number  of  5  j 

i  i 

the  gas  flow  in  the  experiments  was  9.  The  Reynolds  number  was  i  j 

7  j  I 

4.5x10  per  meter.  The  other  operation  made  was  the  straight  j  I 

through  mode.  In  this  mode,  the  storage  chamber  total  pressure  ; 

was  45  bars.  The  total  temperture  was  680°K.  The  Mach  number  of  ,  i 

the  gas  flow  in  the  experiment  was  6.6  and  the  Reynolds  number  was 

1.8x10^  per  meter.  The  duration  of  the  experiments  was  always  4 

milliseconds.  According  to  theoretical  predictions,  at  the  place 

where  the  protuburence  was  installed,  the  boundary  layer  of  the 

flat  plate  would  respectively  have  thicknesses  of  5-3  and  4.4mm.  1 

j 

i  ! 

Measurement  systems.  The  rates  of  surface  heat  flow  for  the  j  | 

flat  plate  were  measured  by  the  use  of  platinum  plated  electrical  |  j 

resistance  temperature  devices.  The  measurement  signal  passed  j 

through  the  thermoelectric  simulation  net  and  was  changed  into  a  j 

signal  which  was  In  direct  proportion  to  the  heat  flow.  After  this  '  i 

was  done  [9],  it  was  recorded  by  means  of  oscilloscope  photography.  ! 

The  data  was  displayed  in  comparison  form.  This  was  done  in  order 
to  get  rid  of  the  error  introduced  by  the  measurement  system  itself. 

That  is  to  say,  before  the  protuburence  was  installed,  experiments 

were  conducted  to  measure  the  flat  plate  surface  heat  flow  rate 

for  the  location  in  question.  After  this  was  done,  a  comparison 

was  carried  out  with  the  surface  heat  flow  rates  measured  after 

the  protuburence  was  installed  and  the  value  of  the  ratio  of  these  j 
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Figure  1.  Model  flat  plate  dimensions  and  symbols, 
two  values  was  solved  for. 

Models .  When  one  compares  convex  protuburences  and  air¬ 
craft  or  spacecraft,  the  dimensions  of  the  former  are  much  smaller. 
If  one  makes  use  of  a  scale  for  geometric  reduction,  the  dimen¬ 
sions  of  a  convex  protuburence  on  a  model  would  be  too  small. 

This  makes  experimental  observations  as  well  as  the  placement  of 
sensors  difficult  to  carry  out.  In  response  to  these  problems, 
we  took  a  flat  plate  on  the  surface  of  the  spacecraft  or  aircraft 
involved  and  in  the  vicinity  of  the  convex  protuburence  Involved 
and  used  this  to  simulate  the  protuburence  on  the  model.  The  Mach 
number  was  also  simulated  with  the  Mach  number  for  the  local  area 
of  the  craft  model  in  which  the  protuburence  was  supposed  to 
appear.  In  this  way,  the  dimensions  of  the  convex  protuburence 
on  the  flat  plate  of  the  model  could  be  much  larger.  The  external 
dimensions  of  the  flat  plate  simulation  used  in  these  experiments 
as  well  as  various  types  of  symbols  are  presented  in  detail  in 
Figure  1.  The  flat  plate  was  installed  with  an  agle  of  attack  of 
zero,  and  in  its  leading  edge  for  20mm  a  2x2mm  horizontal  convex 
trench  was  opened.  This  was  to  facilitate  the  turning  of  boundary 
layers . 

Concerning  the  oil  flow  visualization  technique,  the  process 
of  starting  up  a  shock  wave  wind  tunnel  is  fast.  As  far  as  the 
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reverse  jet  operational  configuration  used  in  these  experiments 
goes,  the  duration  of  the  experiments  was  immediately  stopped  by 
the  sparseness  of  the  reflected  jet  in  the  test  section.  Because 
of  this,  the  process  of  stopping  is  also  very  fast.  Basically, 
surface  oil  flow  flow  spectra  are  formed  by  the  effects  of  gas 
flow  during  the  duration  of  the  experiments.  Although  the  dura¬ 
tion  of  the  experiments  was  very  short,  the  oil  flow  flow  spectrum 
for  the  inside  of  the  area  close  to  the  cylinder  on  which  surface 
there  are  exerted  large  shear  forces  stands  out  very  clearly.  The 
paint  type  material  used  in  the  oil  flow  technique  is  diffusion 
pump  oil  with  a  blackening  powder  added.  Before  the  experiments, 
the  paint  type  material  was  spread  evenly  on  the  surface  of  the 
flat  plate.  After  the  experiments,  the  flow  spectrum  remained. 

Under  conditions  in  which  M  =  4,  the  standard  oil  flow  flow  spec¬ 
trum  is  as  shown  in  photograph  1  (for  the  photographs  in  this 
article  see  Plate  I). 

3.  RESULTS  OF  EXPERIMENTATION  AND  DISCUSSION 

1 .  The  principles  governing  the  distribution  of  heat  transfer 
parameters  in  long  vertical  cylinders.  On  the  basis  of  research 
done  on  the  oil  flow  flow  spectra  of  surfaces  [3]  as  well  as  on 
the  results  of  observations  made  in  these  experiments,  it  is  poss¬ 
ible  to  obtain  the  interference  effects  which  cylindrical  protubur- 
ences  have  with  flat  plate  surfaces.  When  the  length  to  diameter 
ratio  is  L/D  »  5,  it  is  possible  to  consider  the  cylinder  as  un¬ 
limited  in  length.  Concerning  the  long  cylinders  presented  in  this 
article,  their  length  to  diameter  ratios  were  all  equal  to  or 
larger  than  5. 

On  the  center  line  of  flow  over  protuburences ,  there  are  two 
types  of  Mach  number  (M  »  6.6  and  9)  vertical  cylindrical  heat 
transfer  parameter  distribution  curves  which  are  shown  in  Figure  2. 
(.b  is  the  flat  plate  heat  transfer  coefficient  at  measurement 
points  where  there  are  no  protuburences).  The  basic  characteristics 
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of  these  distribution  curves  are  as  follows.  Beginning  from  the 
vicinity  of  the  place  at  which  interference  starts,  the  heat  trans¬ 
fer  parameters  gradually  increase.  After  maintaining  a  steady 
plateau,  these  values  drop.  After  this,  from  their  lowest  point, 
there  is  a  sudden  rise  to  their  peak  values.  After  passing  these 
peak  value  points,  the  values  then  go  back  to  a  monotonously  steady 
decline.  In  order  to  take  the  surface  sections  of  the  neighboring 
columns  and  represent  their  curves  more  clearly,  one  can  use  loga¬ 
rithmic  coordinates,  and  the  area  close  to  the  columns  in  question 
are  expanded  as  shown  in  Pigure  3>  From  this  figure,  it  can  be 
seen  that  the  location  of  the  peak  values  is  at  a  place  from  0.1 
to0'15  from  the  surface  of  the  column  involved. 

The  form  of  the  heat  transfer  coefficient  curves  discussed 
above  can  go  together  with  the  interference  flow  field  structures 
which  can  be  deduced  on  the  basis  of  surface  oil  flow  flow  spectra 
(Plate  1)  and  Schlieren  photography  (photograph  2(a))  In  order  to 
add  to  an  overall  explanation.  There  is  a  group  of  compression 
waves  In  existence  at  the  place  where  the  oblique  shock  waves 
caused  by  the  interference  of  protuburences  and  the  boundary  layers 
of  the  flat  plate  Involved.  The  curves  defined  by  the  coefficients 
of  heat  transfer  begin  their  gradual  rise  at  this  point.  Later, 
these  curves  come  under  the  influence  of  the  oblique  shock  waves 
involved  and,  because  of  this,  they  maintain  a  long  level  config¬ 
uration  in  this  later  section.  Even  later  on,  they  become  subject 
to  the  Influence  of  vortices  at  the  base  of  the  protuburences 
Involved.  Close  by  the  upper  reaches  of  these  vortices,  there  is 
the  transition  area  in  the  direction  of  flow  speed  change.  The 
flow  speeds  are  very  low  and,  because  of  this,  a  trough  section  of 
the  curves  appears.  Close  by  the  base  section  of  the  protuburences 
involved  and  under  the  Influence  of  the  vortices  we  have  mentioned, 
the  coefficients  of  heat  transfer  rise  very  fast.  Due  to  the  fact 
that  the  base  section  of  the  protuburences  involved  are  sharp  right 
angle  turns,  these  vortices  are  not  able  to  directly  exert  their 
influence  on  the  base  section  of  the  protuburences.  Because  of  this, 


7 


m  »  9,  Re  -  4.5x10'/.  meters 


>6.6,  Re  =  1.8x10’/ meters 


Figure  3-  Distribution  of  heat  transfer  coefficients  in  the 
area  near  the  columns  involved  on  the  center  line  of  the  unDer 
reaches  of  round  vertical  columns. 


the  peak  values  are  located  at  a  fixed  distance  from  the  base  sec 
tion  of  them. 


From  the  surface  oil  flow  flow  spectra  (plate  1),  it  is  poss¬ 
ible  to  observe  the  dimensions  of  the  vortices  on  the  surface  of 
flat  plates  involved.  From  the  Schlieren  photography  (Photograph 
2a)',  it  is  possible  to  make  out  the  thickness  of  the  vortices 
involved  in  a  vertical  direction.  The  black  belt  on  the  upper 
boundary  is  a  shear  layer.  This  vortex  is  influenced  by  high 
pressure  coming  from  three  different  directions,  and  this  condition 
is  produced  by  the  downwash  flow  which  adheres  closely  to  the  sur¬ 
face  of  the  columns  involved.  There  is  a  strong  relationship 
between  the  strength  of  the  vortices  concerned  and  the  Mach  number 
of  the  gas  flow.  The  higher  the  Mach  numbers  involved  are  for  the 
gas  flows  concerned,  the  higher  is  the  strength  of  the  vortices 
and,  because  of  this,  the  higher  are  the  peak  values  for  the  heat 
transfer  coefficients. 

If  one  makes  a  comparison  with  two-dimensional  forward  stepped 
heat  transfer  coefficient  distribution  curves  [10],  one  will  find 
similar  configurations.  However,  the  range  of  disturbance  is  much 


smaller  than  that  of  the  two-dimensional  stepped  case  and,  because 
of  this,  the  values  for  the  peak  heat  transfer  coefficients  are 
much  higher. 

Lateral  heat  transfer  coefficient  distribution  curves  for 
vertical  columns  are  shown  in  Figure  4 .  They  are  the  same  as 
those  for  the  upper  center  lines  and  the  peak  value  points  are 
also  located  a  certain  fixed  distance  from  the  surface  of  the 
columns  involved.  This  distance  is  greater  than  the  distance  from 
the  upper  center  line.  However,  the  gradient  in  the  vicinity  of 
the  peak  values  is  smooth  and  gentle.  The  peak  values,  in  this 
case,  are  clearly  lower  than  the  peak  values  for  the  center  line 
of  the  upper  parts.  In  fact,  they  are  almost  only  1/3  of  the 
latter. 

The  horizontal  distribution  curves  for  the  lower  center  of  a 
column  at  a  distance  one  radius  of  the  column  from  it  are  shown 
in  Figure  5.  At  the  intersection  points  on  the  center  line,  the 
coefficients  of  heat  transfer  for  long  columns  and  the  coefficients 
of  heat  transfer  for  non-columner  flat  plates  are  quite  close  to 
each  other.  The  coefficients  of  heat  transfer  toward  the  outside 
decrease  and,  at  the  location  0.2D  from  the  center  line,  they  fall 
to  approximately  half  of  the  flat  plate  values.  Starting  from 
this  value,  they  begin  to  slowly  beat  up  to  their  peak  values. 

After  passing  their  peak  values,  they  very  slowly  drop  down  again. 
These  peak  values  only  have  about  1/7  of  the  peak  values  on  the 
center  line  of  the  upper  parts  of  the  flow;  however,  they  are 
still  several  times  the  values  which  one  gets  for  a  flat  plate. 
Because  of  changes  of  this  sort  are  very  slow,  the  width  of  the 
band  of  heat  increase  and  the  total  amount  of  temperature  increase 
are  still  not  low. 

2 .  An  empirical  formula  for  the  maximum  heat  transfer  coeffi¬ 
cient  value  for  long  vertical  columns.  The  maximum  heat  transfer 
coefficients  for  the  flat  plate  interference  area  appear  on  the 


10 


center  line  of  the  upper  flow.  Cline  [7].  on  the  basis  of  exper¬ 
imental  data,  (2.4  c  M  s  4.44),  and  the  following  empirical  for¬ 
mula  is  obtained 

*"  "  (1) 

In  this  equation,  the  magnitude  of  the  calculated  Reynolds  number 
is  the  immersed  length  of  the  flat  plate.  If  one  makes  a  compar¬ 
ison  of  values  for  the  heat  transfer  coefficients  at  their  highest 
values  calculated  according  to  the  formula  presented  above  in  these 
experiments,  as  far  as  M  *  9  is  concerned,  the  value  is  15.7.  For 
M  =  6.6,  it  is  13.3.  The  deviations  of  46  and '31  from  the  mea¬ 
sured  values  is  quite  considerable.  And  the  deviation  between  the 
compared  values  are  not  the  same  either.  It  appears  that  this 
formula  is  not  suitable  for  high  Mach  numbers. 

Truitt  [8]  takes  the  center  line  flow  and  considers  it  a  two- 
dimensional  flow,  so  that  he  derives  the  relationship  presented 
below, 

WA„  -  ( 1  +  K1M)V5/Re,/’  C2  ) 

On  the  basis  of  experimental  data,  we  get  a  precise  determination 
of  K2  =  25.  A  comparison  of  the  curves  involved  is  shown  in  Fig¬ 
ure  6.  In  the  same  figure,  one  also  sees  the  4B-128  data  from 
reference  [7].  It  can  be  seen  from  this  figure  that  the  Burbank 
[6]  data  and  other  data  of  a  similar  sort  are  wrong  on  the  low 
side.  There  are  two  reasons  for  this.  One  is  the  scarcity  of 
measurement  points  Cround  columns  with  D  =  2.8"  and  an  interval 
between  sensors  of  1").  This  scarcity  makes  it  impossible  for 
accurate  measurements  to  be  made  of  the  actual  peak  values.  The 
second  reason  is  that  because  we  used  a  tin  stainless  steel  shell 
to  measure  heat  flow,  it  was  not  possible  to  correct  for  the  rela¬ 
tively  large  errors  introduced  because  of  lateral  heat  diffusion. 

On  the  basis  of  the  analysis  of  George  [11],  the  error  which  is 
caused  by  laterally  induced  heat  flow  when  one  uses  thin  shells 
for  measurement  purposes  is 
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rate  of  laterally  Induced  heat  flow  _ 
rate  of  surface  heat  flow  a  0 

In  this  equation,  a  is  the  coefficient  of  thermal  conductivity, 
t  is  the  measurement  time.  Q  is  the  rate  of  heat  flow  on  the 
surface.  On  the  basis  of  the  heat  transfer  coefficient  distribu¬ 
tion  curves  obtained  in  these  experiments,  the  place  at  which  one 
gets  the  estimated  peak  value  for  heat  flow,  that  is, 

is  on  the  order  of  10  and  the  value  for  stainless  steel  of  a  = 

2 

0.05  cm  /sec.  The  laterally  induced  heat  flow  error  which  comes 
from  the  use  of  thin  shells  is  capable  of  reaching  a  value  several 
times  the  result.  Due  to  the  fact  that  the  place  at  which  the 
peak  value  occurs,  that  is  ,  is  a  function  of  the  Mach 

number  of  the  gas  flow  involved,  when  this  number  is  reduced,  the 
peak  values  quickly  drop.  Because  of  this,  when  the  Mach  numbers 
involved  are  low,  the  laterally  induced  heat  flow  error  is  small. 

In  Figure  6,  the  error  revealed  in  a  comparison  between  M  =  2.4 
and  2.65  type  numerical  data  and  data  obtained  when  M  *  3-51  and 
4.44  was  small;  this  phenomenon  is  probably  related  concerning 
cases  in  which  one  uses  thin  film  electrical  resistance  temperature 
measuring  devices  to  measure  heat  flows,  due  to  the  fact  that  the 
low  value  Ca^0.008  cm  /sec)  is  extremely  low,  the  reaction  is  fast, 
and  the  measurement  Interval  is  pushed  forward  in  time.  Even  if 
one  takes  his  measurements  at  the  place  at  which  the  heat  transfer 
distribution-  is  extremely  pointed,  the  error  from  lateral  heat 
conduction  will  still  be  very  low. 

3.  The  Influence  of  the  height  of  the  cylinder.  From  Figure 
3,  it  can  be  seen  that,  whether  M  *  9  or  H  =  6.6,  when  the  height- 
to-radius  ratio  is  L/D  >  1,  the  configurations  of  heat  transfer 
coefficient  curves  on  the  center  line  of  the  upper  reaches  of  flow 
for  different  height-to-radius  ratios  are  similar,  and  the  curves 
conform  closely  to  one  another.  When  L/D  <  1,  despite  the  fact 
that  the  configurations  of  the  curves  are  similar,  the  horizontal 
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15*  O  L/0-1  •*{£ 


Figure  4.  Lateral  heat  transfer 
coefficient  distributions  for  long 
colums  (X/D*0  cross  section) 

Key:  1 — top  graph;  2 — meter;  3 — bottom 
graph;  4 — meter;  5 — represents  long  ver¬ 
tical  columns;  6 — represents  long  cylin¬ 
der  swept  back  by  15°;  7 — represents  long 
cylinder  swept  back  by  30°;  8 — represents 
long  cylinder  swept  forward  by  15°;  9 — 
represents  vertical  columns  for  L/D=l 
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dimensions  have  been  foreshortened,  and  there  is  an  obvious  drop 
which  also  occurs  in  the  peak  values  for  the  heat  transfer  coeff¬ 
icients. 


As  was  said  above,  the  heat  transfer  coefficients  for  the  flat 
surfaces  adjacent  to  protuburences  are  strongly  dependent  on  the 
dimensions  and  strengths  of  the  vortices  involved.  When  the 
height  of  a  column  is  smaller  than  its  radius,  the  Y-point  is  not 
far  from  the  apex.  The  high  pressures  in  the  vioinity  of  the  Y- 
point  leaks  off  the  apical  point  early  and  this  causes  a  reduction 
in  both  the  dimensions  and  the  strength  of  the  vortices  involved. 
Because  of  this,  the  differences  between  the  heat  transfer  curves 
for  the  various  height  cases  are  not  great. 

The  curves  which  are  defined  by  the  influence  of  the  height- 
to-radius  ratio  of  cylinders  on  the  maximum  heat  transfer  coeffi¬ 
cients  on  the  center  line  of  the  upper  parts  of  the  flow  are  drawn 
out  in  Figure  7.  The  influence  of  height-to-radius  ratios  indicate 
weakening  and  attenuation  and  inductions  can  be  made  of  the  type 
of  empirical  formula  presented  below,  that  is. 


There  is  a  relationship  between  and  the  Mach  number  of  the 
gas  flow  involved.  In  these  expeirments,  when  M  *  9>  *  1.0, 

and  when  M  =  6.6  then  K.^  *  1.6. 

4 .  The  Influence  of  the  angle  of  Inclination  in  cylinders. 
The  influence  of  sweepback  angle  and  forward  sweep  angles  have  on 
the  curves  defined  by  the  heat  transfer  coefficient  distribution 
along  the  center  lines  of  the  upper  parts  of  flow  are  very  large 
(Figure  8).  First  of  all,  a  sweepback  angle  changes  the  configur¬ 
ation  of  these  distribution  curves.  Even  if  the  sweepback  angle 
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is  15°,  the  distribution  curves  involved  will  already  show  a  change 
from  the  base  seccion  outward  in  the  form  of  a  monotonous  decline; 
however,  there  will  still  be  a  trough.  When  the  angle  of  sweep- 
back  involved  is  30°,  the  trough  will  also  disappear.  After  this, 
angles  of  sweepback  cause  the  areas  influenced  by  interference  to 
be  greatly  reduced  in  size.  The  most  important  thing  is  that  the 
heat  transfer  coefficients  are  greatly  reduced. 

A  sweepback  angle  configuration  is  an  effective  measure  to 
use  in  order  to  improve  the  interference  heating  environment.  It 
is  at  its  most  effective  when  Mach  numbers  are  high.  The  use  of 
configurations  swept  forward  was  found  to  be  without  any  use  at  all 
as  far  as  Improving  the  heating  environment  on  flat  plate  surfaces 
goes.  Due  to  the  strengthening  of  the  high  pressure  downwash, 
there  is  not  only  a  rise  in  the  peak  values  of  heat  transfer  coeff¬ 
icients,  but  also  an  expansion  in  the  range  of  Interference. 

The  influence  which  a  sweepback  angle  configuration  has  on 
the  heating  environment  of  flat  plate  surfaces  can  explain  the 
wave  graphs  which  appear  In  the  Schlieren  photographs  concerned, 
plates  2a,  b  and  c.  If  one  Is  considering  the  case  of  vertical 
columns,  the  bow-shaped  shock  waves  above  Y-points  are  parallel 
to  the  surface  of  the  cylinders,  and  the  gas  flow  behing  these 
waves  are  perpendicular  to  the  axis  line  of  any  round  columns 
that  may  be  Involved.  Concerning  cylinders  with  sweepback  angles 
of  15°  and  30°,  this  same  sort  of  bow-shaped  shock  wave  and  the 
surfaces  of  the  column  make  a  small  included  angle.  (When  the 
angle  of  sweepback  Is  15°,  the  included  angle  is  1°.  When  the  angle 
of  sweepback  is  30° >  the  included  angle  is  2°).  Because  of  this 
fact,  gas  flow  which  follows  the  bow-shaped  shock  waves  is  inclined 
upward.  The  high  pressure  at  the  Y-point  flows  upward,  and  the 
downward  Influence  on  a  flat  plate  is  weakened.  From  the  curves 
In  Figure  8,  it  can  be  seen  that,  when  the  angle  of  sweepback  is 
15°,  the  curves  defined  by  the  heat  transfer  coefficients  involved 
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still  have  troughs.  Because  of  this,  it  is  still  possible  to 
have  a  small,  weak  vortex. 


The  influence  of  inclined  slopes  on  lateral  lines  are  the 
same  as  those  explained  above.  The  only  difference  is  that  the 
influences  are  somewhat  weaker  (Figure  4). 

5.  Gas  flows  and  the  Influence  of  Mach  numbers.  Mach  numbers 

have  an  important  influence  on  the  peak  values  of  heat  flow.  From 

formula  (2)  one  can  see  that  the  influence  of  Mach  numbers  is  in 

1  8 

direct  proportion  with  M  .  As  far  as  lateral  and  lower  flow  heat 
transfer  is  concerned,  the  influence  of  Mach  numbers  have  proper¬ 
ties  which  are  the  same.  Gas  flow  Mach  numbers  also  influence  the 
range  of  areas  of  interference.  When  Mach  numbers  are  high,  the 
range  of  interference  areas  is  large  (Figure  2).  As  far  as  long 
cylinders  are  concerned,  M  ■  6.6,  and  the  area  of  interference  is 
smaller  than  2.5D;  moreover,  when  M  *  9,  interference  areas  are 
larger  than  2.5D.  In  the  case  of  short  cylinders  (L/D  =  0.25), 
these  areas  increase  from  ID  to  1.2D. 

When  gas  flow  Mach  numbers  are  high,  the  weakening  and  atten¬ 
uation  proceeding  down  the  flow  are  relatively  rapid,  that  is  to 
say,  the  influences  of  these  high  Mach  numbers  on  the  lower  parts 
of  the  flow  are  relatively  strong.  Concerning  weakening  and 
attenuation  which  is  caused  by  sweepback  when  gas  flow  Mach  num¬ 
bers  are  high,  the  influence  of  such  swept  back  configurations 
is  strong.  When  M  *  9,  the  peak  value  on  the  center  line  of  the 
upper  parts  of  the  flow  drops  to  one-sixth  of  the  peak  value  for 
vertical  cylinders  when  the  sweepback  angle  is  30°.  Moreover, 
when  M  »  6.6,  it  only  drops  to  approximately  one-fourth. 

4.  CONCLUSIONS 


1.  In  the  course  of  these  experiments,  we  discovered  that 
the  highest  rates  of  surface  heat  flow  within  the  interference 
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areas  of  cylindrical  protuberances  are  not  located  in  the  base 
section  of  the  cylinders  involved  but  in  locations  0.1-0.15D 
from  the  base  sections  of  these  cylinders  on  the  center  line  of 
the  upper  parts  of  the  flow.  Due  to  the  fact  that  the  density  of 
the  measurement  points  and  the  measurement  process  do  not  produce 
lateral  heat  conduction  errors,  the  meausred  peak  values  of  heat 
transfer  coefficient  are  much  higher  than  the  numerical  values 
which  have  been  determined  outside  of  China  from  experiments  and 
calculated  from  formulas.  When  M  =  9,  the  peak  values  for  heat 
transfer  coefficients  reach  values  as  much  as  46  times  higher. 

2.  Swept  back  configurations  are  capable  of  making  obvious 
reductions  in  the  range  of  interference  areas  as  well  as  of  greatly 
lowering  the  maximum  heat  transfer  coefficients.  These  effects  are 
particularly  obvious  when  Mach  numbers  are  high.  Moreover,  forward 
inclined  configurations  are  useless  when  it  comes  to  heating 
environments . 

3.  Gas  flow  Mach  numbers  have  important  influences  on  the 
maximum  heat  transfer  coefficient  values.  These  maximum  heat  trans¬ 
fer  coefficient  values  are  directly  proportional  to  Ml. 8.  When 
Mach  numbers  are  high,  interference  area  ranges  are  expanded  and 
there  is  a  reduction  and  slowing  down  of  the  attenuation  of  the 
lower  part  of  the  flow. 

4.  When  cylinder  heights  are  larger  than  the  corresponding 
cylinder  diameters,  the  influence  of  height  is  extremely  weak. 

When  the  height-to-diameter  ratio  is  small  than  1,  the  influence 
exerted  by  cylinder  height  is  very  large. 

These  experiments  were  completed  in  cooperation  with  the  heat 
transfer  team  and  the  high  pressure  vacuum  measurement  team  of  the 
JP  Laboratory.  The  Schlieren  photography  was  taken  by  comrades 
Liu  Fang  and  Hu  Jin-ming. 
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Abstract 

The  results  of  hent  trnnfer  rate  measurements  in  the  interaction  region  of  a  cylin¬ 
drical  protuberance  is  described.  It  is  found  that  the  maximum  heat  rate  on  the  cen¬ 
terline  does  not  occur  at  the  foot  of  the  cylinder  but  at  a  distance  of  0.1  to  0.15  diameter 
of  cylinder  upstream  from  it.  For  a  cylinder  normal  to  the  surface,  the  maximum  heat 
transfer  rate  is  44i  times  that  of  flat  plate  value  at  Mach  number  9,  and  31  times  at 
Mach  number  6.6. 

The  effect  of  cylinder  sweep  on  both  the  extent  of  interaction  region  and  the  ma¬ 
ximum  heat  transfer  rate  is  significant  particularly  at  high  Mach  numbers.  When  back- 
sweep  angle  varies  from  0s  to  30s,  the  maximum  heat  transfer  rate  decreases  by  a  factor 
of  about  *.  For  forward  sweep,  both  the  extent  of  the  interaction  region  *ud  the  ma¬ 
ximum  heat  transfer  rate  increase. 


THE  PROBLEM  OP  IGNITION  IN  A  LAMINAR  BOUNDARY  LAYER  OVER  A 
HEATED  PLATE* 

Sun  Bing-hua 

(QIng  Dao  Gas  Turbine  Factory) 

Figure  1  shows  a  high  temperature  inertia  heated  plate  for 
a  temperature  T  .  The  combustible  mixture  of  gases  which  flows 
over  this  surface  forms  a  laminar  flow  boundary  layer.  Where 
x  =  xi}  the  combustible  gases 'begin  to  ignite.  This  article  is 
going  to  make  use  of  differential  equations  in  order  to  do  research 
into  this  question. 


Tm.  ft . 


Figure  1.  Ignition  on  a  heated  plate  in  a  laminar  flow 
boundary  layer. 

Key:  1 — speed  boundary  layer;  2 — temperature  boundary  layer; 

3 — heated  plate;  4 — Ignition  point 

BASIC  ASSUMPTIONS 

1.  The  speed  of  the  conflicting  flows  uB  is  very  small. 

2.  The  direct  proportion  between  the  viscosity  u  and  the 
coefficient  of  thermal  conductivity  X  and  the  temperature  T  and 
the  specific  heat  Cp  is  a  constant. 

3.  The  Prandtl  number  Pr  and  the  Schmidt  number  Sc  are  equal 
and  are  constant. 


This  article  was  received  on  April  8,  1978 
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0. 


4.  The  pressure  gradient  ^  = 

5.  The  temperature  of  the  heated  plate  is  TQ  =  cons  t. 

6.  The  speed  of  the  second  degree  two  molecular  reaction 
can  be  represented  as  the  reaction  speed 

Wi  -  K*pViexP  (“  jf). 

In  this  equation  Kq  is  the  speed  constant  for  the  reaction, 
f  is  the  corresponding  concentration.  R  is  the  normal  gas  con¬ 
stant.  E  is  the  reactivity,  i  represents  the  i-th  type  of 
component . 

BASIC  EQUATIONS 

If  we  make  use  of  the  Howarth-Dorodnitsyn  transformation,  it 
is  possible  to  take  the  basic  equation  for  a  compressible  air 
flow  and  transform  it  into  the  basic  equation  for  incompressible 
air  flow,  that  is 


dx  dy 
du  du 

ut^  +  vw~ 


9T  .  ar 

■I  -  +  V - - 

6x  dy 

df .  dh 

*  d7  +  ‘'  d7 


a** 

V~W 

v,  a»r  +  (-jf) 

Pr  dy2  C, 

Vm  d^fi  (  E  \ 

S"dp  ”  Kut>^  “P(~  RTJ 


Cl) 


In  these  quations,  Qx  is  the  heat  value.  The  corresponding  bound¬ 
ary  conditions  are 
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y_o,  «-  *-o,  r-  r„  l^-o, 

ay 

(2) 

y  — ►  co  ,  «—■««,  T  —  T« ,  /,  —  /,» 

From  (1)  and  (2),  if  we  carry  out  an  integration,  it  is  easy  to 
obtain  the  differential  for  the  energy  equation  and  the  diffu¬ 
sion  equation  involved.  This  takes  the  form  given  below,  that  is 


L* -  /,*,  -  g  (f ). + k.  J>«*  (-  i) « 


(3) 

C4) 


In  order  to  make  the  induction  easy,  take  (3),  (4)  and  convert 
them  to  a  dimensionless  form.  In  this  way,  make 


V  - 


u 


F 


fit  /i- 


c- 


r  * 


In  these  equations,  respectively  represent  speed,  temp¬ 

erature  and  the  thickness  of  the  concentration  boundary  layer, 

In  this  way,  (3),  (.4)  are  transformed  into  the  form 


dx 


r.  (  dd  \ 
Pr  umfl  a?  /• 


PQiiKtiPt1]ttxV  (  ) 

cr(r,  —  r.)«. 


C5) 


V.  I 

1  dF\ 

i  a. 

(-nr) 

ScUmT 

vac  /' 

•  (/.-  —  /,o)«- 

(6) 


In  these  equations,  In  represents  a  differential,  that  is 


l‘-(J 

LYexpf- 

(— 1  )| 

(.7) 

,0  RTt 
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APPROXIMATION  CROSS  SECTIONS 


On  the  basis  of  a  Poblhausen  selection,  as  far  as  a  flat 
dP 

plate  defined  by  -3—  -  0  is  concerned,  the  speed  cross  section  is 


U  —  —  —  2?;  —  2 t;1  4-^',  “  4" 

Urn  O 


(8) 


On  the  basis  of  the  boundary  conditions  below,  we  select  the 
temperature  cross  section,  that  is 


y  —  0,  7  —  To, 


(  E  \ 

Qyp  *  Pf  A  n».oo/,o®^P  ^  'r  T 


3yJ 

ar  ajr 


VmCf 


y  -  P,  T  -  Tm,  dy  —  dyt  ~  0 


We  can  introduce  a  non-dimensional  temperature  deformation  para¬ 
meter,  that  is 

t  P„AV«/;<>e*P  I  -  ,77-) 

-  -  '  n~  •  3* 

df(  l't  —  Tm  )fm 


II  CO 


(.9) 


If  we  make  use  of  a  cross  section  determining  precise  values 

for  tempeature  and  consisting  of  more  than  four  quantities,  the 

result  is  g  _  (  _  i.  ( 12  — //)f  -  —  f' +  —  (4  + //)f5 

T,  —  Tm  6  2  2 


-1(6+  //)?*,  B  -  4 

6  p 


(10) 


The  ignition  condition  is  x  -  x,  and  ( f— )  —0 

1  •  oy  ® 

From  (1G)  it  is  possible  to  obtain  the  temperature  deformation 
parameter  for  the  moment  of  ignition  which  is 

H(x,)  -  12.  (ID 
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"The  significance  of  the  temperature  deformation  parameters 
is  that,  when  H  increases  with  x,  the  temperature  fields  in¬ 
volved  give  rise  to  violent  deformations.  Obviously,  this  is 
a  result  due  to  the  liberation  of  heat  in  chemical  reactions. 

The  temperature  fields  involved  are  drawn  out  in  Fig.  2. 


The  selection  of  concentration  profiles  by  the  use  of 
similar  methods  is  as  presented  below: 


F  -  - k  -  1  —  -1  (12  —  f)C 

/-  -  /«  6 

2  2 

6  T 


(12) 


In  this  equation,  c  is  the  concentration  deformation 
parameter,  that  is 


(13) 


1.0 


1  -  H  =»  0,  2-H-6,  4  -  H  -  12  4  -  W  =.  18. 


If  one  is  speaking  from  the  point  of  view  of  an  inertia 
heated  plate,  the  boundary  condition  for  concentration  is  ^ _0 
At  this  time,  on  the  basis  of  the  concentration  cross 
seciton,  there  must  be 

COO  »  12.  (14) 

In  this  type  of  situation,  f  is  a  function  of  x,  or  it  is  a 
function  of  H.  On  the  basis  of  the  definition  for  H,  Equation 
(9)  as  well  as  (.14)  show  the  relationship 
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« 1  —  HM  (/I  —  1), 
M  —  C,Tm 


(15) 


In  this  case,  it  is  necessary  to  point  out  that  the  concen¬ 
tration  boundary  which  is  used  by  Dooley  is  fiQ  »  const.  Morever, 
that  used  by  Toong  is  —  o  •  The  solutions  for  these  are 

both  numerical  y  solutions,  Obviously,  as  far  as 

inertia  heated  plates  are  concerned,  the  second  type  of  boundary 
condition  is  relatively  reasonable. 

APPROXIMATE  SOLUTION  EQUATIONS 
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As  far  as  the  boundary  condition 
is  a  function  of  x.  One  can  take 


(*L)  - 

'  dy  /• 


is  concerned, 
the  temperature 


deformation  parameter  and  change  it  so  that  it  is  written  in  the 

form  .  E  . 

Pr  Q.iKufxJ', mexp[— 

C,(Tt—Tm)vm  (16) 


H  -  nNY,  n 


If  one  takes  the  corresponding  cross  section  and  substitutes  it 
into  Equation  (5),  then  one  can  obtain 


-j-  IMi  +  A2H)  1  -  -  ~~  (12  -  H)  +  -2=_«JW  *  /. 

dx  r  PtUmfi  Pttlm 


g 

In  this  equation  and  h2  are  both  functions  of  r  = 
when  they  are  differentiated,  one  gets  the  following 


when  r  $  1 


when  r  >  1 
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1080r 


(.17) 

and 

(18) 

(18a) 


These  calculations  demonstrate  that,  within  the  range  r  =  0.7  -  1.8 
the  difference  between  the  two  equations  can  be  completely  ignored 
sind  not  taken  Into  consideration  for  purposes  of  these  calculations. 
In  this  way,  one  can  make  use  of  (18)  later  because  of  the  fact 
that  it  is  relatively  easy  to  use. 


The  differential  I 


becomes 


(JL9) 


If  one  takes  (17)  and  carries  Its  simplification  a  step  further, 
one  obtains 
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The  non-dimensional  coordinates  in  these  equations  are 


X 


(20c) 


Equation  (20)  is  nothing  else  than  an  approximate  solution  equa¬ 
tion  for  H  and  r.  It  is  still  necessary  to  supplement  the  situa¬ 
tion  with  one  more  equation.  Prom  the  definition  of  H,  one  can 
get  the  equation  (1 6),  that  is 

H  —  nX1?  —  r.X-r‘81 


On  the  basis  of  the  Karman-Pohlhausen  approximate  solution,  the 
thickness  of  the  speed  boundary  layer  is 


If  one  substitutes  the  equation  above,  then  he  gets 


//  -  m  Pr ;[  i  -  n\t (a  —  i ) iv. 


(.21) 


This  is  nothing  else  than  the  supplementary  equation  which  we 
required.  The  initial  conditions  for  this  are  as  follows: 

When  i-0,  W-0,  r-r,--4=.  C22) 

</Pr 

The  condition  which  we  are  considering,  that  is, 
obtained  in  the  following  way.  We  took  Equation 
and  substituted  into  Equation  (20),  that  is 


r»  ”  T— 
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and  the  denominator  i~*o 


In  the  equation  above, 
the  numerator  must  — 


when  x  — o,  H-*o 


o,  and  in  such  a  case, 

-i  —  34  Pc  r,j  -  o 

A, 


one  gets 


This  equation  is  an  equation  which  solves  for  the  ratio  r  of  the 
speed  boundary  layer  thickness  and  the  temperature  boundary  layer 
thickness  when  no  combustion  reactions  are  taking  place.  The 
value  of  this  is  r,  —  _ 1 

VpT 

Q 

It  should  be  explained  that  it  is  only  when  one  selects  —  =1 
that  it  is  possible  to  make  calculations  of  (19)  with  relative 
ease.  The  analysis  involved  is  the  same  type  as  that  used  in  the 
cases  of  (18)  and  (l8a).  The  error  which  is  given  rise  to  by 
the  use  of  this  method  is  very  small.  In  this  way,  we  can  then 
forget  about  Equation  (6).  Otherwise,  it  would  be  necessary  to 
take  Equation  (6)  into  consideration,  and  this  would  make  the 
problems  involved  much  more  complicated. 
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Abstract 


Thd  problrau  of  ignition  and  combustion  in  laminar  boundary  layer  orer  a  heated 
plate  are  studied  by  means  of  an  approximate  method  analogous  to  that  developed  by 
Karm&n  and  Pohlhausen  in  their  investigation  of  the  separation  of  laminar  boundary 
layer  ever  a  plate. 

The  basic  partial  differential  equations  are  transformed  into  an  ordinary  differen¬ 
tial  equation  of  the  first  order  for  the  temperature  deformation  factor  B,  and  the  cal¬ 
culation  of  integral  in  is  a  prerequisite  to  resolving  the  equation.  The  introduction  of 
factor  B  is  covenient  in  determining  the  ignition  distance  since  B  equals  12  at  the 
ignition  point  on  the  plate 
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